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The  subject  of  differential  equation  constitutes  a  large  and 
very  important  branch  of  modern  mathematics.  From  the  early  days 
of  Calculus  the  subject  has  been  an  area  of  great  theoretical 
research  and  prat i cal  applications,  and  it  continues  to  be  so  to 
dale.  ?hr-rorore,  a  natural  question  arises:  How  does  one  obtain 
useful  information  from  a  differential  equation.  The  answer  is 
essentially  that  if  it  is  possible  to  do  so,  one  solves  the  differ¬ 
ential  equation  to  obtain  a  solution,  if  it  is  not  possible,  one 
uses  the  theory  of  differential  equations  to  obtain  information 
about  the  solution.  In  the  following  text  we  will  investigate 
some  ordinary  differential  equations  in  self-adjoint  form  and 
when  possible  determine  the  behavior  of  solutions  without  actually 
solving  the  differential  equation. 
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The  following  definitions  will  be  usedi 

SELF-ADJOINT  EQUATION  is  a  second  order  linear  differential 

equation  of  the  form  [/(x^']1  +  p(x)y  -  0  ,  where  r(x) 
and  p(x)  are  continuous  and  r(x)  >  0  on  an  interval  a  £•  x  tb. 

A  BOUNDARY- VALUE  PROBLEM  is  a  problem  which  consist  ofi 

1.  a  second  order  homogenuous  linear  defferential  equation 

of  the  form  jr(x)y'T  +  jp(x)  +  \q(x)]  y  =  0,  where  p,  q,  (1) 
and  r  are  real  valued  functions  such  that  q  has  a  continuous 
derivative,  p  and  r  are  continuous,  q(x)>  0  and  r(x)>  0  for 
all  real  x  on  an  interval  a^-x^b,  and  X  is  a  parameter 
independent  of  x, 

2.  two  supplementary  conditions 

A1y(a)  +  A2y’(a)  =  0 

Biy(b)  +  B2y'(b)  =  0 

where  A^ ,  A2,  B^ ,  and  B2  are  real  constants  such  that  A ^  and 
A2  are  not  both  zero  and  Bj  and  B2  are  not  both  zero. 

A  STURM-LIOUVILLE  PROBLEM  consists  of  a  differential  equation  of 
the  form  [jr(x)y'3'  +  jp(x)+  q(x)]y  =  0,  and  boundary  condi¬ 
tions 

Ajy(a)  +  B1y* (a)  =  0 
A2y(a)  +  B2y*(a)  =  0. 

A  SINGULAR  STURM-LIOUVILLE  EQUATION  has  any  one  of  the  following 
conditions « 

1.  the  interval  is  finite  and  p(x)  vanishes  at  one  or  both 
endpoints, 

2,  the  interval  is  infinite  or  semi-infinite. 
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A  REGULAR  STURI'-LIOUVILLE  EQUATION  has  an  interval  which  is  closed 


and  finite,  that  is,  for  every  real  x  the  at’Xib. 

EIGENVALUES  are  the  values  of  \  in  the  differential  equation 
{r(x)y,3'  +  [p(x)  +  q(x)]y  =  0,  for  which  there  exist  non¬ 
trivial  solutions  of  the  problem. 

EIGENFUNCTIONS  are  the  corresponding  nontrivial  solutions  them¬ 
selves. 

The  unique  graph  of  the  solutions  to  an  initial  value  problem  in 
an  appropriate  region  is  often  referred  to  as  a  solution 
path. 

The  isolated  critial  point  (0,0)  is  called  a  saddle  point  if 

there  exists  a  neighborhood  of  (0,0)  in  which  the  following 
conditions  holdj 

1.  There  exist  two  paths  which  approach  and  enter  (0,0) 
from  opposite  directions  as  t-»a>  and  there  exist  two  other 
paths  which  approach  and  enter  (0,0)  from  different  opposite 
directions  as  t-*-oo  . 


2.  In  each  of  the  four  domains  between  any  two  of  the  four 
paths  in  (1)  there  are  infinitely  many  paths  which  are  arbi¬ 
trarily  closed  to  (0,0)  but  which  tend  away  from  (0,0)  as 
t-^  +oo  and  as  t  ->  -oo  . 


In  the  study  of  self-adjoint  equations  we  need  a  well-known 

theorem  on  point  sets  which  is  the  "olzano-’.Viestrass  Theorem. 

Now  supoose  E  is  a  set  of  ooints  on  the  x-;ix3s.  A  ooint  x  is 

*  o 

called  a  limit  point  of  E  if  there  exist  a  sequence  of  distinct 

po'nts  x'  ,  x„,  x?,  ...  of  S  such  that  lim  x  -  x  .  The  Tolzano- 
J  n  -*cd  n  0 

V/iostrass  Theorem  states  that  every  infinite  set  E  has  at  least 

one  limit  point. 

All  differential  equations  can  be  put  in  self-adjoint  form. 

In  the  discussion  that  follows  we  will  be  investigating  the  zeros 
of  the  solutions  of  self-adjoint  equations. 

Theorem  1.  Let  f(x)  be  a  solution  of  (r(x)y']+  p(x)  =  0 
whose  first  derivative  f'(x)  exists  and  is  defined  on  aixsb. 

If  f(x)  has  an  infinite  number  of  zeros  on  a±x  f  b  then  f(x)  =  0 
for  all  x  on  a£Xib. 

Proof:  Since  f(x)  has  an  infinite  number  of  zeros  on  ja,^ 
then  by  the  Bolzano-Wiestrass  Theorem  the  set  of  zeros  has  a  limit 
point  xQ,  where  xQ  6  [a,bl.  Thus  there  exist  a  sequence  ^x^  of 
zeros  which  converge  to  xQ  (with  xn  /  xQ).  Since  f(x)  is  con¬ 
tinuous  then  lim  f(x),  where  x-*x  through  any  sequence  of  points 
x  -*a> 

on  (a, 6/  .  Let  x  ~*x0  through  the  sequence  of  zeros  .  Then 

lim  f(x)  =  0  =  f(x_).  Since  f'(x)  exists,  we  have 

x-4xn  0 

o 

f'(x  )  =  lira  f(x)  -  f(a) 

X  ~>x  x  -  x^ 

0  o 

where  x-*xQ  through  the  sequence  jx^  .  Now  for  such  points 

f(x)  -  f(xQ)  we  have  f*(x)  =  0.  Therefore  by  the  Uniqueness 
x  "  xo 
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theorem  for  Differential  Equation?*  f(x)  is  a  solution  of  a 
self-adjoint  equation  such  that  f(xQ)  =  f'(x  )  =  0  for  all  x 
on  a  f-  x  c-  b .  Q .  5 .  D . 

AB5L*5  FCRITJLA.  Given  two  solutions  f(x)  and  g(x)  of  a 
self-adjoint  equation  on  an  interval  aix'-b,  then  gor  all  x  on 
a-Jx-b  we  have  r(x)  [f(x)g'  (x)  -  f'(x)g(x)J  =  K. 

Proof:  Since  f(x)  and  g(x)  are  solutions  of  a  self-adjoint 

equation  on  a  ix  i  b,  we  have 

£r(x)f*(x|  +  p(x)f(x)  =  0  (2) 

and  jr(x)g'  (x)j'  +  p(x)g(x)  =  0  .  (3) 

Multiplying  equation  (2)  by  -g{x)  and  equation  (3)  by  -f(x)  and 
adding  the  resultant  equations  we  obtain 

f(x)  (r (x)g'(xjl  ~g(x)  [r(x)f,(xj]1  =  0. 

Integrating  from  a  to  x  and  using  integration  by  parts  we  have 

*  r* 

f(t)r(t)g'(t)  -  / r(t)g'(t)f'(t)dt  -  g(t)r(t)f *(t) 

+  f r{x)f '(t)g*(t)dt  =  0 

-'a 

which  redueces  to 

r(x)  [f(x)g'(x)  -  f '  (x  )g(  x  )j  =  r(a)  jf  (a)g'  (a)  -  f*(a)g(a)J 
The  right  hand  member  is  a  constant  and  by  letting  K  represent 
this  constant  we  have  r(x)  f(x)g'(x)  -  f'(x)g(x)  =  K  Q.E.D. 

Lemma  1«  If  two  solutions  f(x)  and  g(x)  of  a  self-adjoint 
equation  have  a  common  xero  then  they  are  linearly  dependent. 
Conversely,  if  f(x)  and  g(x)  are  linearly  dependent  solutions, 
neither  identically  zero,  then  if  one  of  them  vanishes  at  x  =  xQ 
so  does  the  other. 
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Proof:  For  the  first  part  we  shall  apply  Abel's  Formula. 

Let  xq  t  |a»b)  be  a  common  zero  of  f(x)  and  g(x).  Letting  x  -  xq 
we  obtain  k  =  Oj  hence,  r(x)  jf(x)g*(x)  -  f'(x)g(x)j  -  0  for  all 
x,-  (a,b].  Since  we  have  assumed  r(x)>  0  on  a  ^  x  -t  b  this  means 
that  the  quantity  in  the  brackets  must  be  zero.  Put  the  quantity 
in  the  brackets  is  \V(f*g)(x)  which  is  the  V/ronskian  for  f  and  g. 
By  '.Vronskian  definition  of  linearly  dependent,  if  the  determinant 


YJ  Jf(x) ,  ,g(xj]  = 


=  0 


f(x)  g(x) 
f'(x)  g'(x) 

then  the  solutions  f(x)  and  g(x)  are  linearly  dependent. 

Second  part.  Since  f(x)  and  g(x)  are  linearly  dependent  on 
airx-b,  then  there  exist  constants  and  c2,,not  both  zero, 
such  that  c-^fCx)  +  c^gCx)  =  0.  By  our  hypothesis,  neither  f(x) 
or  g(x)  is  zero.  Hence  if  Cj  =  0  then  c^g(x)  =  0,  but  g(x)  /  0 
which  means  c2  must  be  zero.  But  this  is  a  contradiction  to 
c1  and  c2  both  being  zero.  Hence  ^  0  and  c2  /  0.  By  letting 
x  =  xQ  then  if  f(xQ)  =  0  then  g(xQ)  =  0.  Q.E.D. 

Theorem  2.  (Sturm  Separation  Theorem).  Let  f(x)  and  g(x) 
be  linearly  independent  solutions  of  [r(x)y'J  +  p(x)y  =  0  on  the 
interval  a-xib,  then  between  any  two  consecutive  zeros  of 
f(x)  there  is  precisely  one  zero  of  g(x). 

Proof:  Let  x^  and  x2  be  two  consecutive  aeros  of  f(x)  on 
(a,bj.  Then  by  lemma  1  we  know  that  g(x^)  /  0  and  g(x2)  /  0. 
Assume  that  g(x)  has  no  zero  on  x^x^x^  Since  the  solutions 
f(x)  and  g(x)  have  continuous  derivatives  on  a,b  then  the 
quotient  f(x)/g(x)  has  a  continuous  derivative  on  x^x^Xg. 
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Furthermore  f (x)/g(x)  -  0  at  the  endpoints  of  this  interval. 


’’once  by  Rolle's  ^heorem  there  exist  ^  where  ^  c.  such  that 


lut  jf (x)/g(x )]  -  and  since  f(x)  and  r ( x )  are 

L  J  [g(x)]  2 

linearly  independent  on  b,  we  have  jf(x)/g(xjj  /  0 

on  x^<x<x^.  But  this  is  a  contradiction.  Hence  g(x)  has  at 
least  one  zero  on  x.  4  x<x^. 

How  suppose  g  has  more  than  one  zero  in  x^  px<x2.  Let  x^ 
and  x^  be  the  two  consecutive  zeros  of  g.  By  interchanging  f  and 
g  in  the  above  discussion,  this  shows  that  f  must  have  at  least 
one  zero,  say  x^,  in  the  open  interval  x^4x^<x^.  Then  we  would 
have  x1<^x^x2  meaning  that  x1  and  x2  are  not  consecutive  zeros 
of  f.  This  is  a  contradiction  to  our  assumption  that  x^^  and  x2 
were  consecutive  zeros.  Therefore  g  has  precisely  one  zero  in 
the  interval  x1<x<x2<  Q.E.D. 

Theorem  3.  (Sturm  Comparison  Theorm)  Let  f(x)  be  a  solution 
of  /r(x)y^j  +  p(x)y  =  0  and  g(x)  be  a  solution  of  r(x)y'  +  j^(x)  =  0. 
Let  r(x)>0  and  have  a  continuous  derivative  and  let  p(x)  and 
Pj ( x )  be  continuous  and  p1(x)>p(x).  If  Xj  and  x2  are  consecutive 
zeros  of  f(x)  on  [a,b] ,  then  g(x)  has  at  least  one  zero  at  some 
point  on  the  interval  x^4x<x2 

Proofs  Assume  g(x)  does  not  have  a  zero  on  the  interval 
x^t  x4x2<  By  our  hypothesis  we  have 

[r(x)y'J>  +  p(x)y  =  0 
{rCxJy'J1  +  p1(x)y=  0 
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for  all  x  e  |~a,b].  Now  we  multiply  the  first  equation  by  g  ( x )  and 
the  vi.  1  equation  by  f(x)  and  subtract.  The  resul tant  equation  i 
)  i  r(x  )  f '  (  x  )]'  -  f(x)  [r(x)g*  (x)J  =  [p^x)  -  p(x)j  f (x)g(x) 

"he  left  hand  side  of  the  equation  yields 

•-(  x  )  |r(x)  f  '(:■:)]  -  f(x)  jr(x)r,’(x)]  -  jr{  x  )  [f *  ( x  )r  ( x )  -  f(x)g'(xjj. 

There  fore 

(r(x)  Q '(x)g(x)  -  f(x)g'(x)]]'  =  Ppx(x)  -  p(x)]f(x)  g(x). 

Integrating  from  x.  to  x0  we  have 
jf{r(x)  [f  ’  (x)g(x)  -  f(x)g’  (xjj)'  dx  =J[P1(x)  -  p(x)]  f(x)g(x)dx 
or 

r(x)  [f '  (x)g(x)  -  f(x)g'(xj]  =  ^[PiCx)  -  p(x)]  f(x)g(x)dx 
Since  f(x1)  =  f(x2)  =  0  then  the  above  equation  becomes 
r(x2)f '(x2)g(x2)  -  r(x1)f’(x1)g(x1)  =  ^0^! ( x )  -  p(x)J  f(x)g(x)dx 
3y  our  hypothesis  r(x2)>0.  Now  since  f(x2)  =  0  and  f(x)>0  on 
x1^x<x2,  we  have  f'(x2)c0.  Also  since  g(x)>0  then  g(x2)^  0. 
Hence  r(x2)f ' (x2)g(x2) £  0.  In  a  similiar  manner,  we  have 
r(x^ ) f ' (x1 )g(x^ ) >  0 .  Therefore  the  left  hand  side  of  equation 
(4)  is  not  positive.  Also,  by  our  hypothesis  p^(x)  >p(x)  implies 
that  p1(x)  -  p(x ) >  0 .  So  the  right  hand  side  of  equation  (4)  is 
positive.  This  is  a  contradiction  to  the  assunption  that  g(x) 
does  not  have  a  zero  on  the  interval  x^Z  X4  x2>  Q.E.D. 

Theorem  4.  (Osgood  Theorem)  Suppose  p(x) >  0  and  p'(x)>  0 
on  the  interval  0  4x<r  rnand  let  y(x)  be  an  arbitrary  solution  of 
y’’  +  p(x)y  =  0.  If  x=a  and  x=b  are  two  consecutive  zeros  of 
y '  (x )  then  |y(b)|  i  j y ( a )|. 


(4) 
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Proof : 


Suppose  p’(x)-0  on  the-  interval  aSrx^b,  then  p(x) 
is  consent  -nd  Jy(b)j  -  jy(a)|  .  Tf  p'(x)/0  and  y(x)  la  a  -ela¬ 
tion  of  y''  +  p(x)  -  0  then  the  equation  becomes 

y'  '  (x)  +  p(x)y(x)  -  0  (f) 

'■'ultinlyi  ng  equation  (f)  by  2y'(x)  and  i  n  terrat  i  ng  a  ho  resultant 
equation  over  the  interval  a,b  v:e  obtain 


y'2(x)  +  j  v(x)  [y2(x)J  '  cx  =  0 

O  ^ 


"he  first  term  on  the  left  is  sero,  threfore  our  equation  becomes 

b  9 

x)  y-'(x)  '  dx  =  0. 


Using  integration  by  parts  we  have 


fp(x)  y2(x)  •  dx  =  p(b)y2(b)  -  p(a)y2(a)  -  /p' (x)y2(x )dx  =  0 
J0  . 

p' (x)y2(x)dx 


"o 

or 


(6) 


Suppose  y2(b)>y2(a)f  then  p(b)y2(b)  -  p(a)y2(b) 

p(b)y2{b)  -  p(a)y2(a), 

2  2 

then  p(b)y(b)  -  p(a)y  (b)  is  greater  than  both  the  left  and  right 
members  of  equation  (6).  It  follows  that 

P(b)y2(b)  -  p(a)y2(b)  4  p(b)y2(b)  -  p(a)y2(a) 

~  p(a)y2(b)  7  -  p( a)y2( a) 
p(a)y2(b)  4  p(a)y2(a) 

But  this  contradicts  our  assumption  y  (b)4  y  (a).  Therefore 
y2(h).^y2(a)  and  |y  (b)  |  *  j  y  ( a)  j  on  the  interval  06.x-£oo.  Q.E.D. 

Theorem  j.  Let  r(x)  be  positive  and  suppose  that  r(x)  and 
p(x)  are  continuous  on  the  interval  04x400.  If 

I  1  dx  =  +  od  ,  and  (o) 

J  r(x) 
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(3) 


non  ov- 


>(x)dx  -  +  oo 


-ry  solution  y(x)  of  equation  (1)  vanishes  infinitely 


often  on  'one  interval  U  x  c<x>  .  Siiniliarly,  in  the  interval: 


/  rT’xT  dx  =  +  ®  •  and 

c 

j  p(x)dx  =  +  oo 


then  every  solution  y(x)  of  equation  (1)  vanisnes  infinitely 
often  on  the  interval  O^x^l. 

Theorem  6.  Every  nontrivial  solution  of  /r(x)y^J*  +  p(x)=0 
has  at  most  a  finite  number  of  zeros  on  the  interval  x  <  qd  ,  if 


dx  4  +  cd  ,  and' 


f&T  * 

/  p(x)dx  4  M 
Ja 


where  M  is  any  positive  constant. 

Theorem  7-  (Lienard-Levinson-Smith  Theorem)  A  differential 

p 

equation  d  x  +  f(x)dx  +  g(x)  =  0 
dt^  dt 

(K 

where  f,  g,  and  F  defined  by  F(x)  =  I  f(u)du  and  G  defined  by 
vt  J° 

G(x)  =  j g(u)du  are  real  functions  having  the  following  properties t 

1.  f  is  even  and  continuous  for  all  x. 

2.  there  exists  a  number  x  V  0  such  that  F(x).4  0  for  0^x<.x_, 

o  o 

F(x)>0  and  monotonic  increasing  for  x>xQ.  Furthermore 
v/e  have  F(x)->ao  as  x-*co  . 

3.  g  is  odd,  has  a  derivative  for  all  x,  and  g  is  such  that 
g(x)  >  0  for  x  >0. 
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4.  !(>:)  -*no  at'  x  '>;D  . 

’’vjn  the  differential  equation  d'x  +  f(x)r]x  +  g  ( x )  =  0 

dx  -  dt 

has  a  unique  nontrivial  periodic  soultion.  In  other  words  the 

equivalent  autonomous  system  cte  =  y  ,  =  -  f(x)y  -  g(x) 

dt  dt 

has  a  unique  closed  path  in  the  xy  plane. 

Froofs  Recall  an  autonomous  system  is  a  system  in  which  the 
independent  variable  appears  only  in  the  differentials  on  the  left 
members  and  not  in  the  functions  on  the  right.  In  the  above  system, 
t  is  the  independent  variable.  Considering  the  differential  eoua- 

/)  p 

tion  d * x  +  f (x )dx  +  g(x)  =  0  where  f(x)=E(x  -  1),  E  constant,  and 
dt?  dt 

g(x)~x  we  have  the  well  known  Van  der  Pol  equation 

d2x  +  S(x2-l)dx  +  x  =  0. 
dt^  dt 

It  follows  that  F(x)  =  f f(u)du 

=  /*E(u2-l)du 

=  E(u£  -  u)  x 
3  0 

=  E(x£  -  x) 

3 

and  G(x)  =  ^  g(u)du 

=  J*1  udu  =  x£  . 

Checking  to  see  if  our  properties  hold  we  have 

1.  since  f(-x)  =E(x2-l)  =  f(x),  then  the  function  f  is 
even  and  clearly  f  is  continuous  for  all  x. 
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2.  F(x)  -  F(x  /3  -  x)  =  x2/3(x2-x).  Here 

F('J'T)  =  0 

F ( x  )  <0  if  x<  \J  3^ 

F(x)>  0  if  x  >  vT* 

F(x)  is  nonotone  increasing  and  F  >oo 

as  t-iao  . 

3-  Since  g(-x)  =  -x  =  -g(x),  it  follows  that  the  function 
g(x)  is  odd.  Also  since  g*(x)=l,  then  the  derivative  of  g 
is  continuous  for  all  x.  It  is  obvjously  true  that  g(x)">  0 
for  x>  0. 

2  2 

Therefore  the  Van  der  Pol  ecuation  d  y  +  E(x  -l)dx  +  x  =  0 

dt2  dt 

satisfies  all  the  properties  of  the  Lienard-Levinson-Smith  Theorem 

and  the  equation  has  a  unique  nontrivial  periodic  solution. 

Furthermore  the  equivalent  autonomous  system 

dx  =  y  d^  =  E(  1  — x  )y  -  x 

dt  dt 

has  a  unique  closed  path.  Reducing  the  above  autonomous  system 

to  a  differential  equation  of  the -paths,  we  have, 

dy  dx  =  E ( 1 -x  )y  -  x  y 
dt  dt 

Using  the  chain  rule  we  have 

d 2L  =  E(l-x2)  -  x  .  Q.E.D. 

dx  y 

Now  we  will  discuss  how  to  put  differential  equations  in 
self-adjoint  form.  Given  a  linear  differential  equation  of  the 
second  order 

a(x)y"  +  b(x)y*  +  c(x)y  =  0 


(13) 
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whore  a(x)">0  and  a(x),  b(x),-and  c(x)  are  continuous  on  the  inter¬ 
val  fa,b],  this  equation  can  be  put  in  self-adjoint  form  by  the 
follow! ng  procedure;  First  multiply  both  members  of  the  equation 
by  the  function 

1  e 

a(T7 


liext  set  r(x)  = 


'm 


dx 


ffiBH3dx 


and  p(x)  =  c(x)  e 
a(x) 

The  self-adjoint  equation  will  be  l[r(x}y'3  +p(x)y=0.  A  necessary 
and  sufficient  condition  that  equation  (13)  be  self-adjoint  is  that 
[a(xj]=  b(x)  on  a<  x  £  b. 

Using  the  above  procedure,  we  shall  put  the  Bessel,  Hermite, 

and  Laguerre  equations  in  self-adjoint  form.  Consider  the  Bessel 
2  2  2 

equation  x  y"  +  xy'  +  (x  -n  )y  =  0,  where  n  is  a  constant.  It 
follows  that 


b  ( x  )1  dx 
la(x). 


a(x) 


=  1  e 

2 


1  dx 
x 


=  1 
x 


Multiplying  both  members  of  the  Bessel  equation  by  1/x,  we  have 

xy"  +  y*  +  [&2-n2.)/x]y  =  0, 

2  2 

Therefore  r(x)=x  and  p(x)=x  -n  ,  and  the  self-adjoint  form  of 

x 


the  Bessel  eauation  is  xy'  '  +  (x-n_)y  =  0 

x 


is  a 


Consider  the  Hermite  equation  y 
constant,  we  have  J^xj J  dx 

alx) 


'  -  2xy  *  +  2ny  =  0  where  n 

-*/xdx  2 

e  J  =  e"x 


(14) 
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Mu 1 M plying  each  term  by  e_x  we  find  r(x)~e~x  and  p(x)-?ne  x  . 

"herefore  the  self-adjoint  form  of  the  Hermit  equation  is 

r_  2  2 
je'X  y'J  '  +  2ne“X  y  =  0  . 

Consider  the  Laguerre  cauation  xy'*  +  (l-x)y'  +  ny  =  0, 

where  n  is  a  constant.  We  have 

\  Mx )  dx  f(  1-x )  dx 

J a(x)  J  x 

1  e  =  1  e  =  e_x  . 

a(x )  x 

_x 

Multiplying  each  term  by  e  we  obtain 

e“xxy * '  =  (l-x)e-xy*  +  ne-xy  =  0  . 
Therefore  the  self-adjoint  form  of  the  Laguerre . equation  is 

|xe-xy^|  '  +  ne  xy  =  0  . 

Now  we  are  going  to  give  examples  of  how  to  investigate 
oscillation  of  solutions  of  differential  equations  without 
actually  solving  the  differential  equation. 

The  Sturm  Separation  Theorem  essentially  says  that  the  zeros 
of  one  of  two  linearly  independent  solutions  of  a  self-adjoint 
equation  separate  the  zeros  of  the  other  solution.  That  is,  the 
number  of  zeros  on  an  interval  of  a  self-adjoint  equation  is  about 
the  same  as  the  number  of  zeros  of  the  other  solution.  For 
example,  we  can  show  that  between  any  two  consecutive  zeros  of 
sin  x  +  ^cos  x  there  is  one  zero  of  sin  x  -  cos  x.  We  know 
by  the  Sturm  Separation  Theorem  that  the  zeros  of  one  linearly 
independent  solution  of  a  self-adjoint  equation  separates  the  zeros 
of  the  solutions  of  the  other  self-adjoint  equation.  First  we 
shall  construct  our  equations  in  self-adjoint  form  and  apply 


(15) 
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Sturm's  Separation  Theorem.  In  the  first  function  we  let 

v  -  sin  x  +  \(?cos  x  ,  then 

v *  *  —  —  sin  x  -\2cos  x,  therefore 

v"  +  v  =  0  is  our  self-adjoint  equation  with  r(x),  p(x)  =  1. 

Tn  the  second  function  we  let 

z  -  sin  x  -  cos  x  ,  then 

z"  =  -  sin  x  +  cos  x  . 

Here  we  have  another  self-adjoint  equation  z1'  +  z  =  0  .  Thus 
v(x)  and  z(x)  are  two  linear  independent  solutions  of  y"  +  y  =  0. 
3y  applying  Sturm's  Separation  Theorem  we  know  between  any  two 
consecutive  zeros  of  sin  x  +  \)5cos  x  there  is  one  zero  of 
sin  x  -  cos  x. 

Given  two  self-adjoint  equations,  what  can  one  say  about  their 
solutions?  We  find  that  the  Surm  Comparison  Theorem  can  best 
answer  this  question.  The  theorem  compares  the  rates  of  oscil¬ 
lation  of  soultions  of  the  two  equations. 

We  now  consider  the  two  differential  equations 

2 

x  y*'  +  xy'  +  y  =  0  and  y*'  +  y  =  0  and  investigate  which  possesses 
the  more  rapidly  oscillating  solution  on  the  interval  (1,  oo) . 

We  know  by  Theorem  3  that  if  two  equations  are  in  self-adjoint 
form  and  Pj(x)>p(x)  then  the  solutions  of  the  self-adjoint  equa¬ 
tion  whcih  contains  Pj(x)  will  oscillate  more  rapidly.  The  second 
equation  is  already  in  self-adjoint  so  all  that  remains  is  to  put 
the  first  equation  in  self-adjoint  form  and  apply  Theorem  3. 

Putting  the  first  equation  in  self-adjoint  we  [xy]j  '  +  ly  =  0. 

x 

Hence  our  two  equations  become 
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ixy 


J 


'  +  ly  -  0  and  y '  '  +  y  =  0.  Since  <^1  for  all  x  in  1  <  x  <  ao  , 
x  x 


then  by  the  Iturra  Comparison  theorem  the  solutions  of  y"  +  y  =  0 

will  oscillate  more  rapidly  than  the  solutions  of  ^xy'J  •  +  ly  =  0 

x 

Suppose  we  are  given  two  differential  equations 

x^y ' '  +  xy*  +(x^-l)y  =  0  ,  and 
xz"  +  z '  +  xz  =  0 


(16) 

(17) 


with  nontrivial  solutions  u(x)  and  z(x)  respectively,  and  both 
equations  vanish  at  x=l.  Now  we  determine  which  equation  will 
vanish  first  after  x=l.  This  problem  can  best  be  solved  by  first 
putting  the  two  equations  in  self-adjoint  form,  and  applying  Theo¬ 
rem  3*  We  know  by  Theorem  3  if  u(x)  and  z(x)  are  solutions  of 
two  self-adjoint  equations  then  the  differential  equation  with  the 
larger  second  term  will  oscillate  more  frequently.  Our  objective, 
therefore,  would  be  to  reconstruct  the  two  equations  and  compare 
the  coefficient  of  the  second  terra  of  each  equation.  Both  equations 
vanish  at  x=l,  so  we  must  find  which  equation  oscillates  with 
the  highest  frequency  beginning  after  x=l.  The  respective  self- 
adjoint  forms  of  equations  (16)  and  (17)  are 

rfj  '  +  (x-l/x)y  =  0  ,  and 


z '  *  +  lz  =  0  . 

X 

Considering  only  the  second  term  of  both  equations  for  x>l  we 

have  1  ^  x-1^  .  Therefore  by  Sturm's  Comparison  Theorem  between  any 
x  x 

tv/o  consecutive  zeros  of  z"  +  (l/x)z  =  0  there  is  one  zero  of 

/ly']  '  +  (x-l)y  =  0.  Since  solutiions  of  equations  ( 1 6)  and  (17) 
Lx  -*  x 
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both  have  a  common  "oro  at  x=l,  then  between  x=l  and  the  next 
zero  of  z'*  +  (l/x)z  =  0  there  is  one  zero  of  equation  (16). 
Therefore  solutions  of  equation  ( 1 6 )  oscillates  with  higher  fro- 
quency  than  solutions  of  equation  (17),  and  xy*'  +  xy'  +(x  -)  =  0 
will  vanish  first  after  x=l. 

In  our  discussion  thus  far  we  have  considered  the  case  where 

p(x)>0  and  have  seen  that  a  nontrivial  solution  of  jr(x)y'J  '+p(x)y 
can  have  infinitely  many  zeros  on  the  interval  ^a,b].  One  may  ask 
the  question,  what  can  be  said  about  the  nontrivial  solutions 
if  p(x)^0.  In  the  following  discussion  we  shall  see  that  if 
p(x)^r  0  then  no  nontrivial  solution  of  ^(xjy’j'  +  p(x)y  =  0  can 
have  more  than  one  zero  on  fa,b]  . 

Let  u(x)  be  a  nontrivial  solution  of  r(x)y*  +  p(x)y  =  0, 
and  let  xQ  be  a  zero  of  u(x)  so  that  u(xQ)=0.  Since  u(x)  is  a 
a  nontrivial  solution,  then  u’(x)/o,  therefore  we  have  two  cases, 
u'(xo)>0  or  u'(xo)<0. 

Case  1,  Suppose  u'(xo)>0,  then  u(x)  is  positive  over 
some  interval  to  the  right  of  xQ.  By  the  hypothesis,  the  equation 
is  (r(x)yiJ  '  +  p(x)y  =  0  ,  or 

(r(x)y •  (x >7  *  +  p(x)y(x)  =  0  . 

Using  variable  substitution  we  have 

^(xJu’CxX)  '  +  p(x)u(x)  =  0  ,  or 
i(x)u'{xj]  ’  =  -p(x)u(x)  . 

We  are  given  p(x)^  0.  If  p(x)=0  then  \r(x)u' (x)]  =0  and  we  have 
the  desired  zero.  If  p(x)<0  then  {r(x)u'(x)]  is  a  positive 
function  on  the  same  interval.  This  means  that  the  slope  u’(x)  is 
increasing  so  u(x)  cannot  have  a  zero  to  the  right  of  xQ.  By 
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analogous  reasoning,  u(x)  cannot  have  a  zero  to  the  left  of  xQ. 

Case  2.  Suppose  u'(x)<0  then  r(x)u'(x)  is  a  negative 

function.  This  moans  that  the  slope  u*  (x)  is  a  decreasing  function 

and  u(x)  cannot  have  a  zero  to  the  left  of  xQ.  Similarly  u(x) 

cannot  have  a  zero  to  the  right  of  xQ. 

Therefore,  if  p(x)tO  then  a  nontrivial  solution  of 

(r(x)y'_j'  +  p(x)y  =  0  can  have  no  more  than  one  zero  on  Sa,bJ. 

Now  consider  the  differential  equation  y'*  +  p(x)y  =  0, 

where  p  is  contunuous  and  p(x)'>0  on  a5rX-£b.  Let  pm  denote  the 

2  2 

minimum  value  of  p(x).  We  will  show  that  if  p  \  k  H  then 

^  0 

(b-a)2 

every  solution  of  y"  +  p(x)y=0  has  at  least  k  zeros  on  aSxib. 

First  we  consider  the  two  equations 

y"  +  p(x)y  =  0  ,  and  (18) 

y"  +  k2fr2  y  =  0  .  (19) 

(b-a)2 

We  will  examine  the  solutions  of  equation  (19)  and  apply 

Theorem  3-  A  solution  of  equation  (18)  is  sin  x  and  a  solution 

of  eauation  (19)  is  sin  (k _ \  where  k  is  a  multiple  of  b-a.  By 

\b-a; 

2  2 

assumption  p(x)>k _  ,  so  by  Theorem  3  the  solutions  of  equa- 

(b-a)2 

tion  (18)  oscillates  v/ith  higher  frequency  than  the  solutions  of 

equation  (19).  The  solution  of  equation  (19)  is  sin  k  If  x  where 

b-a 

2  2 

k  is  a  multiple  of  b-a.  Therefore  y"  +  k  ffr  y  -  0  has  at  least 

(b-a)2 

2  2 

k  zeros.  Since  p(x) >  k  n  then  by  Tehorera  3,  between  any 

(b-a)2 

two  zeros  of  equation  (19)  there  is  at  least  one  zero  of  equation 


1? 


(18).  Since  equation  (19)  has  at  least  k  zeros  therefore 
y'*  +  p(x)y  =  0  contains  at  least  k  zeros. 

The  Sturm  Comparison  Theorem  is  very  helpful  in  providing 
information  about  equations  of  the  form  y’’  +  p(x)y  -  0.  Con¬ 
sider  the  equation  y'*  +  p(x)y  =  0  where  p  is  continuous  on 
aixtb  such  that  Ookp(x)<,  M.  Let  u(x)  be  a  solution  of 
y''  +  p(x)y  =  0  having  consecutive  zeros  at  and  x2 
(where  /X^<  b).  Then  f?  /  -  x,  /  ip  . 

"  1  ^  2  1  jir 

Consider  the  following  equation  y''  +  My  =  0  where  v  is 
a  solution  such  that  v(x1  )=0.  Then  y=  I  sin  (x-x.)^]?  .  Since 

1  /vr  1 

M>p(x)  then  by  Theorem  3  there  will  be  at  least  one  zero  of 
equation  y' '  +  lly=0  between  any  two  consecutive  zeros  of 
y*'  +  p(x)y=0t  that  is,  there  will  be  a  zero  of  y"  +  My=0 
between  any  two  zeros  and  x2>  Since  x^  is  a  zero  of  y’’+p(x)y=0 
then  the  zero  of  y"  +  My=0  will  be  at  x,+  7?  and  x0>x1+JT  , 

1  W1*  2  1  TIT' 


hence  x0-x,  >  IT 

2  1' fa¬ 


llow  consider  y* '  +  my=0  where  w  is  a  solution 


such  that  w(x2)=0.  Then  y=  J  sin(x2~x)ffir.  Since  p(x)?m,  then 

T™T" 

by  Theorem  3  there  is  at  least  one  zero  of  y’’  +  p(x)y=0  between 
any  two  consecutive  zeros  of  y"  +  my=0.  By  assumption  Xj  and  x2 
are  two  consecutive  zeros  of  y"  +  p(x)y=0,  therefore  a  zero  of 
y'*  +  p(x)y=0  will  be  at  x.  and  a  zero  will  be  at  X.+  ft  where 

1  1  (W 


Xj+_3L>  Hence  x2<xx+  TT  ,  and  x^-x.  <£  TT  ,  s0  that 

{W1  ^ST  4  sfrU 
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4V'X2“>:1  <  • 

^  1  Nln? 

Applications  of  the  '•  lurm  Comparison  Theorem  can  al.-o  bo  made 
to  the  "enrol  equation  for  specific  values  of  the  constant. 

For  example  let  us  examine  the  Becnel  equation  for  oscillation  of 
its  solutions.  Let  y  (x)  be  a  nontrivial  solution  of  the  fennel 

J  p' 

equation  on  the  positive  axis.  If  0  q  p  x  |  then  every  interval 
of  length  rr  contains  at  least  one  zero  of  yp(x).  If  p=|-  then  the 
distance  between  successive  zeros  of  yp(x)  is  exactly  fr  and  if 
p  then  every  interval  of  length^?  contains  at  most  one  zero 
of  yp(x) . 

By  Theorem  3,  if  y(x)  and  z(x)  are  two  nontrivial  solutions 
of  y"  +  p1(x)y=0  and  y'  '  +  p2(x)y=0  respectively,  and  if 
P1(x)>P2(x),  then  between  two  successive  zeros  of  z(x)  there  is 
at  least  one  zero  of  y(x).  Furthermore,  given  a  differential 
equation  y* '  +  p(x)y'  +  q(x)y=0  it  can  be  put  in  self-adjoint 
form  z'*  +  q1(x)z=0  just  by  a  simple  change  of  the  dependent 
variable.  One  refers  to  the  equation  y' '  +  p(x)y'  +  q(x)y=0  as 
the  standard  form  of  the  equation,  and  z''  +  q^(x)z  as  the  normal 
form. 

In  solving  this  problem  we  take  the  Bessel  equation  in  its 
standard  form  and  put  it  in  its  normal  form.  Then  take  the  dif¬ 
ferent  values  of  p,  substitute  them  into  the  equation  in  its  normal 
form,  and  apply  Theorem  3  by  comparing  the  coefficient  of  the 
second  term.. 

For  equation  z”  +  z=0  a  solution  u(x)  is  sin  x.  Since  sin  x 
has  vanishing  points  at  sin  nb>  where  n=0,  1,2,  .  .  •,  then 


-  *r; 
Sp.V/’- 


between  any  two  success ive  r,oros  of  sin  x,  the  length  of  the 

interval  is  it  . 

Consider  the  Ccssel  equation 

x2y’*  +  xy’  +  (x2-p2)y=0  (20) 

where  p  is  a  parameter.  3y  making  the  transformation 

y u(x)x-;2  ,  this  reduces  equation  (20)  to 

u"  +  jjL  +  (l-4p2)ju  =  0  .  (21) 

4x2 

Case  1.  0^04  If  p=0  then  eq^uation  (21)  becomes 

u*  •  +  1*1  +  1_  lu  =  0  (22) 

4x2J 

Therefore  if  yp(x)  is  a  solution  of  equation  (21)  then  yp(x) 

is  a  solution  of  u"  +  Jl  +  1 _ lu  =  0.  Let  z(x)  be  a  solution  of 

4x2-> 

z  *  *  +  z  =0.  The  coefficient  of u  in  equation  (22)  is  greater  than 
the  coefficient  of  z  in  z*’  +  z  =0.  Therfore  by  Theorem  3  between 
any  two  successive  zeros  of  z"  +  z  =0  there  is  at  least  one 
zero  of  a  solution  of  equation  (22).  We  previously  noted  that 
any  two  successive  zeros  of  z(x)  are  at  intervals  'H' distance 
apart,  therefore  in  any  interval  of  length  ^  there  is  at  least 
one  zero  of  yp(x).  Now  consider  the  condition  0<;p^f.  Using 
the  coefficient  of  u  in  equation  (21)  and  substituting  the  value 
of  the  parameter  p  we  have 

1  +  ( l-4p2)  >  1  +  l-4(£)2 

4x2  4x2 

i-4p2  ?  i=m 

4x2  4x2 
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>  o 

4X2 

fhe  can  f  fi  ciont  of  u  in  eouation  (21)  is  1  +  ( l-4o")  ,  whore 

4x2 

O 

( l-4p  ~ )  0,  and  the  coefficient  of  z  in  s'1  +  z  --0  is  1.  I' 


follows  that  1  +  ( 1  -4p ~ )  y  1,  a: •/!  by  Theorem  3  between  any  two 

4x2 

successive  zeros  of  a  solution  of  z''+z=0  there  is  at  least  one 

zero  of  equation  (21).  But  the  zeros  in  z''+z=0  ?.re  fh  distance 

apart;  therefore,  for  0  ^  p  f  in  equation  (21)  there  is  at  least 

one  zero  in  any  interval  of  length/]1. 

Case  2.  For  p=|  we  have 

u' '  +  fl  +  l-4p2  u  =  0 
4x2  J 


u"  +  Fl  +  1-0  1  u  =  0 
L  4x2J 


u*  •  +  u  =  0. 

The  successive  zeros  of  u"  +  u=0  occur  at  the  same  intervals  as 

the  successive  zeros  of  z"  +  z  =0.  Therefore,  when  p=|  the 

successive  zeros  of  equation  (21)  occur  at  intervals  of  length  IT. 

Case  3.  Considering  only  the  coefficient  of  u  in 

equation  (21)  we  have  1  +  l-4p2  /  1  +  l-4(?)2 

o  o 


4x‘ 

l-4p2  /  1-4(4) 
4x2  4x2 

1-4.P2  o  . 

4x2 


4x' 

2 


Here  the  coefficient  of  u  is  negative.  In  an  earlier  example 
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1 


wo  tv'.o'.vd  'hat  if  n(x)  <  0  In  equation  (1)  i.hen  no  nontrivial 
novation  o-n  have  '.ore  tlvin  one  zero  on  the  interval  a  t  x  i-b. 


1  r/lrt.n; 

V  I,"  2 


>1  for  p>’  ,  then  by  "nooren  3  bo  tween  any  two 


err.  >  vo  roa  Oi 


+  x  0  there  is  at  least  one  '.ero  of  ooua- 


t  i  on  (t'l).  Tut  for  p  >  ;• ,  aqua  ion  (:li)  has  at  most  one  aero  on 
any  interval.  The  successive  zeros  of  z’’  +  z  =0  are  Tt  distance 
apart .  Therefore  for  p  > 1  on  any  interval  TT  distance  apart,  there 

1I-.  r.  +-  ..  ^  ’ll  f  I  1  I  L  t  I  I'  \  ill  A 


is  at  most  one  zero 


of  u"  +  ^1 


u  =  0  . 


From  our  applications  of  Theorem  3  (Sturm  Comparison  Theorem) 
we  know  if  we  are  given  a  differential  equation 

y"  +  p(x)y  =  o  (2 

where  p(x)  is  positive  on  the  interval  0  £  x  -c  oo  ,  then 
p(x)>p(0)>0  and  every  solution  of  equation  (23)  vanishes  infinitely 
often  on  the  interval.  What  can  one  say  about  the  amplitudes  of 
the  oscillation.  Do  they  increase,  decrease,  or  remain  the  same 
as  x  increases  on  the  interval.  W.  Osgood  investigated  this 
problem  and  developed  the  Osgood  Theorem  (Theorem  4).  He  dis¬ 
covered  that  under  the  conditions  stated  above  the  amplitudes  of 
the  oscillations  of  solutions  of  equation  (23)  never  increase  as 
x  increase  on  the  interval  0^tx<od  . 

2  2  2 

Now  let  us  test  the  Bessel  equation  x  y'’+xy'+(x  -n  )y=0 
for  oscillations  on  the  interval  1  <  x  <  ao  .  We  know  from  our 
previous  discussion  that  the  self-adjoint  form  of  this  equation  is 


We  also  know  by  Theorem 


[ xy'J  '  +  (x-njy  =  0 

if  (ldx  =  +  cd  ,  then  every  solution  y(x) 

J  X 
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.  Ml  3*] 


vanishes  infinitely  often  on  t. 


•  »  • 

!,::e  interval  1  <  x  <.  <d  .  ;u t  j  ldx 


does  not  exist  and  the  lost  fails.  We  ai no  k now  from  our  urevi o 
investigation  of  the  "'enrol  equation  that  for  an  interval  of 
length  If  and  certain  values  of  n  we  do  have  oscillations,  J'ak- 
ing  the  traon forma. tin  y -x  'u,  the  Bessel  equation  becomes 

u"  +  fl-4n2-l  lu  =  0 


"low  divide  the  r.-nris  into  intervals  of  length  ft1  .  If  0  4  n  <1 , 
then  every  interval  has  at  least  one  zero  of  a  particular  solution 
y(x).  Therefore  every  solution  y(x)  of  the  Bessel  equation  in 
the  normal  form 

u"  +  I! -4n2-l  u  =  0  , 

-  4x2  -I 

v/ith  0<  n<  vanishes  infinitely  often  on  the  interval  1<  x<  to 
with  the  length  of  each  sub-interval  equal  to  TK 

Our  work  on  oscillations  has  prepared  us  for  a  brief  study 
of  another  areas  The  Sturm-Liouville  Problems.  Our  first  concern 
is  a  study  of  a  special  type  of  two  point  boundary-value  problem. 
There  are  two  important  cases  of  the  Sturm-1 iouville  problem, 
those  in  which  the  supplementary  conditions  are  either  of  the  form 

y(a)=0  ,  y(b)=0  ,  or 

y’(a)=0  ,  y'(4)=0  . 

Consider  the  following  boundary-value  problems 

(xy')’  +  ( 2x2  +  >x3)y  =  0  , 


with  conditions 


3y(l)  +  4y* (1)=0  ,  and 
5y(D  -  3y'(i)=o  . 


Equation  (25)  is  of  the  form  [r(x)y2]  '  +  £p(x)+q(xj]y  =  0 
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2  2 

where  r(x)-x,  p(x)=2x^,  and  q(x)--\x'J;  and  the  conditions  are  of 
the  form  A-^yCa)  +  Bjy’Ca)  =  0  ,  and 

A2y(h)  +  32y’(b)  -  0  . 

no  solve  this  problem  v:e  must  find  a  function  f  which  satisfies 
’noth  the  d  i  fferen  t  i  al  equation  and  the  supplementary  conditions. 

One  solution  to  any  problem  of  this  type  is  the  trivial  solution 
w  such  that  w(x)~0  for  all  values  of  x.  Even  though  this  is  a 
solution  to  our  problem,  one  must  asmit  this  trivial  solution  is 
not  very  useful.  Therefore  we  focus  our  attention  on  the  search 
for  nontrivial  solutions.  During  our  investigation  we  shall  see 
that  the  existence  of  such  nontrivial  solutions  depends  upon  the 
value  of  the  parameter 'K  in  (r(x)y^  ’  +  [p(x)+>q(x)J  y  =  0.  To 
illustrate  this,  let  us  consider  the  Sturm-Liouville  problem 

y*  *  +  y  =  0  (26) 

y(0)=0  ,  y(7?)=0  . 

Vie  have  three  cases:  X=0  ,  >  <  0  ,  or  >  >  0.  In  each  case  we  shall 
first  find  the  general  solution  of  the  differential  equation  and 
then  attempt  to  determine  the  two  arbitrary  constants  in  the 
solution  so  that  the  supplementary  conditions  will  also  be  satis¬ 
fied. 

Case  1.  V=0.  For  this  case  the  differential  equation 

y  *  *+  y=0  reduces  to  y"=0.  Hence,  the  general  solution  is 
y=c^+c2x.  Applying  the  first  supplementary  conditions  y(0)=0, 
we  have  c^+CgTT  =0.  Since  c1=0  then  c2  must  be  zero.  Therefore 
in  order  for  our  solution  to  satisfy  our  conditions,  we  must  have 
c^=c2=0.  Our  solution  then  becomes  y(x)=0  for  all  values  of  x. 

Therefore  for  the  parameter  )^=0,  the  only  solution  of  y"+Ay=0 
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is  !,he  trivial  solution. 

Case  2.  V<,0.  The  characteristic  equation  of  equation  (26) 

P  p  V 

is  n '  +  A  -  0.  It  follows  that  m"  =  -r 

m  =  +  ^ 

oir.ce  V,0,  the  roots  are  real  and  unequal.  Let  a  =  NT-  A  , 

then  the  general  solution  of  equation  (26)  is 

ax  ,  „  -ax 
y  -  c^e  +  c2e 

Applying  the  first  condition  y(0)=0  we  have,  c^+c2=0.  By  the 
second  condition  y(T))~0  we  have,  c^ea^  +  c2e  a^  =0.  In  order  for 
the  solution  y=c^eax+c2e~ax  to  satisfy  the  conditions  y(0)=0  and 
y(n)=0,  the  constants  c^  and  c2  must  satisfy  the  system  of  equations 
c1+c2=0  and  c1ea"+c2e“a^  =0.  Obviously  Cj=c2=0  is  a  solution  of 
this  system  but  this  is  only  the  trivial  solution.  Now  we  find 
the  non  zero  values  of  c1  and  c2<  Recall  that  a  necessary  and 
sufficient  condition  that  n  homogenuous  linear  algebraic  equations 
in  n  unknown  have  a  nontrivial  solution  is  that  the  determinants 


For  the  equality  to  hold,  (a)  must  equal  zero.  That  is,  for  the 
nontrivial  function  y=c^eai^+c2e~arr  to  satisfy  the  conditions 
y(0)=0  and  y  (ir)=0,  we  must  have  a=0.  But  a=’f-A  which  implies 
that  X  =0.  By  assumption  A^O,  therefore  for  this  case  there  are 
no  nontrivial  solutions. 
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Case  3«  V>0.  .  ’  ace  >>0,  then  the  roots  of  the 

characteristic  equation  are  the  conjugate  complex  numbers  +  ■( f i < 
'therefore  the  rcneral  solution  of  equation  ( 2o )  is 

y  =  c^sin  \~j?x  +  CgCos  >7*x. 

Applying  the  first  condition  y(0)=0,  v/e  ha-'o  c^sin  O+c^cos  0=0. 
Since  cos  0=1  then  c^=0.  Applying  the  second  condition  y(/f )  -0 , 
we  have  c^sin  fjF  4-  c2cos  \'f^=0.  Since  c^-0  the  above  equation 
reduces  immediately  to  c^sin  f>tr=0.  fhis  equation  can  be  solved  in 
either  of  two  ways;  set  c^=0  or  set  sin  \]7f  =0.  However,  if  c^=0, 
we  have  already  established  that  c2=0,  then  the  solution  reduces 
to  the  trivial  solution,  c^  must  not  equal  zero.  So  we  must  have 
sin  'f>rt'=0.  But  sin  KtT  =0  only  if  k  is  a  positive  integer.  In 
order  for  equation  (26)  to  have  a  nontrivial  solution  of  the  form 

y=c1sin  '&X+  c^cos 

satisfying  the  conditions  y(0)=0  and  yCffJ^O  we  must  have  )k=n2, 

where  n=l,  2,  3,  ...  .  In  summary,  if  >t:0  then  the  Sturm- 

Liouville  problem  consisting  of  y,,+  V=0  and  conditions  y(0)=0, 

y(fT)=0  does  not  have  a  nontrivial  solution.  However,  if  V?0 

then  a  nontrivial  solution  can  exist  only  if  X  is  one  ot  the 

values  given  by  \=n2  where  n=l ,  2,  3*  •••  •  Furthermore,  from 

the  general  solution  y=c^sin  ^x+c2cos  M"px  we  find  the  non- 

trivial  solutions  corresponding  to  X=n  are  given  by 

y  =  c  sin  nx  where  c  is  an  arbitrary  nonzero  constant.  In  other 
n  n 

words,  the  functions  defined  by  c^sin  x,  c2sin  2x,  c^sin  3 x,  . .., 
where  c^,  c2,  c^,  ...  ,  are  arbitrary  non  zero  constants,  are 
nontrivial  solutions  of  the  problem. 
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A  :> ’  urm-Li ouv ille  equation  can  h  re  an  interval  v;hi  eh  is 
finite,  r^r.i-inf  ini  te ,  or  infinite.  By  the  definitions  at  the 
bo.-- inn  in."  of  the  text,  the  equation  is  classified  either  as  a 
regular  or  singular  vturm-Liouville  equation.  Two  immediate 
examples  of  singular  S  lurm-T.-iouvil  le  equations  are  the  Laxuerre 
and  fennel  eou.a  ions,  how  consider  the  bom idary-value  problem 
yf,  +  >y=0,  with  supplementary  conditions  y(0)=0  and  y ( tT )  =  0 . 

^his  enuation  is  a  S turm-Liouville  equation  with  r(x)=l,  q(x)=l, 
and  p(x)=0.  The  interval  is  finite  and  the  functions  r(x)  and 
a(x)  are  positive  and  continuous,  therefore  the  equation  is  a 
regular  Sturm-Liouville  equation. 

The  investigation  of  y' '+Xy=0  and  supplementary  conditions 
reveal  that  the  values  of  X  for  which  there  exist  nontrivial 
solutions  are  X=n  where  n=l,2,3»  •••  •  By  definition  6,  these 

values  are  are  called  eigenvalues  of  the  problem.  By  definition 
7,  the  eigenfunction  of  the  problem  are  the  corresponding  non¬ 
trivial  solutions  y=cnsin  nx  where  n=l,2,3t  . ..,  and  the  cn  are 
arbitrary  nonzero  constants.  Eigenvalues  and  eigenfunctions  play 
an  improtant  role  'n  many  equations*  One  such  equation  is  the 
Schroedinger  wave  equation  for  the  harmonic  oscillator  in  quan¬ 
tum  mechanics.  This  equation  is 

d2w  +  8TT2m(E-^kx2)w  =  0 

dx2  h2 

where  E  is  the  total  energy,  h  is  a  constant,  and  the  solutions 
w(x)  are  known  as  Schroeding  wave  functions.  Using  the  equation 
k=4n  rav  to  eliminate  the  constant,  the  equation  becomes 


(2?) 
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d?W  +  Pifx  (e-?;>2W)w  -  0  . 

dx2  ■  h2 


,;’ho  ohyrically  admi  ns  ini  n  :olu  1 5  on:;  of  '.hi::  -oca  t  ’  on  are  those 


satisfying  the 

cond 

it  ions  w  ->uo 

as  ) 

X|  >00 

end 

jl 

|2dx 

=  1. 

oc 

,0hcse  solutions 

of 

the  Schroodi 

nror 

wave 

hnc’.i(.,r.s  are  cal'1  0  a 

the  eigenfunctions  and  exist  only  when  E  has  certain  special  values 
called  eigenvalues.  If  we  change  the  independent  variable  to 


1 — 

u  =  2fT  ym  x  ,  or 
x 

_  =  x2  , 

kr2vm 


then  equation  (28)  becomes 

d2w  4  2vm  +  8  2m  ( E-vu2h\ w  =  0  ,  and 
du2h  h2  '  2  ' 

d2w  +  2j  E-u2vh\w  =  0  , 
du2  2  ' 

d2w  +  (2E  -  u2]w  =  0  . 

du2  lvh  ' 

Equation  (29)  is  a  Sturm-Liouville  equation  with  \=2E  -  u2. 

vh 

The  solutions  of  equation  (29)  satisfy  w-*0  as  (u|->  co  and 

(  I  wl2du  =  2\]  vnP  . 

J.  od  h 

Equation  (29),  except  for  the  notation  has  exactly  the  form  of 

another  Sturm-Liouville  equation 

d2z  +  ( 2p  +  1  -  x2)z  =  0, 
dx2 

where  p  is  a  constant.  Physicists  are  only  interested  in  solutions 


=  0  ,  and 
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of  this  equation  that  approach  zero  as  Jx  |  ^ro  ,  therefore  we  shall 

only  consider  those  solutions.  Now  we  take  equation  (30)  and 

transform  it  into  a  well  known  equation,  the  dermite  eouation. 

,2 

"0  S’mplity  equation  (30)  a  new  independent  variable  t=ye-;X 

is  introduced.  This  transforms  equation  (29)  into  the  eouation 

y ' ' -2xy '+2py=0  which  is  the  dermite  equation  with  the  self-adjoint 

_}-2  .  _x2 

form  discussed  earlier  as  e  "  y'*  +  \e  y  =  0  .  Recall  the 
dernite  eouation  has  a  two  form  recursion  formula 

a„ ,  „=  -2(p-n)  a.. 

^+2  Tn+T7Tn+2l  n 

and  the  formula  generates  two  independent  series  solutions: 

y.  (x)  =  1  -  ox2  +  22p(p-2)x^~  -  ...  ,  and 

2!  4! 

y2(x)  =  x  -  2jp_-ljx3  +  23(p-1)(p-3)x5  -  ... 

3S  51 

-l-x2 

Mow  y1(x)e  ’  — >0  as  (xj-^co  if  and  only  if  the  series  y^(x) 
terminates,  that  is,  if  and  only  if  the  parameter  has  one  of  the 
values  0,2,4, ... 

Proof:  IF  PART.  Using  L'Hospital's  rule  we  have 

lim  d/dx  y.(x) 

x  -xp  -3 - 

SL  o 

dx  ei-x* 

=  lim  ( -2px+224p(p-2)x3  -  236p(p-2 ) (p-4)x3  +...)  1 

x-^co  4!  6!  X  2 

xe2X 

=  lim  (-2px  +  22p(p-2)x3  -  23p(p-2) ( p-4)x3  +...)  — — 9  . 

X-^CD  3!  5 •’  xe*x 


Since  the  denominator  goes  to  infinity  faster  than  the  numerator, 
then  the  equation  goes  to  aero  as  ( xf — * gd  . 


ONLY  IF  FAR?. 


"’his  portion  of  the  proof  will  be  shown  using 


the  concept  of  contranosi  tive .  Assume  p  /  0,2,4,...  and  show 
does  not  approach  cero.  y.  (x)  has  the  form 

~  o””*  "*■ 


y1(x) 

an+2~ 


2n 


2  (_p-n)  _  __ 
Tn+TTTn+2 ) 


with  its 


^n  ,  and 


coefficients  determined  by 

the  condition  aQ=l.  Let  ei,X  =J^b2 nx2n 


where  b2n=  1  ,  therefore 

2nn! 

yi ( x )  =  a2nx2n  =  aQ+a2x2+a4x4+  ...  +  a2nx2n+  **' 

-  ~  otr  p  E  2n  * 

e>x~  °2nx  W  +V  +  •••  +  b2nx  +  •" 

To  show  y^(x)  does  not  approach  zero  as  x  go  ,  it  would  be 


sufficient  to  show  a2n>  b2n  if  n  is  large  enough.  Observe  that 


a2n+2  =  -2  (-p-2n) _  ,  and 

a2n  (2n+l)(2n+27 

h2n+2  =  1 

bn  2(n+l7 


Therefore 

a2n+2x'a2n  _  -2(p-2n)  .  2(n+l) 

b2n+2/bn  '  (2n+l)(2n+2> 


For  n  sufficiently  large  we  have,  -2( p-2n)  — *2.  This  implies  that 

2n+l 


a2n+2  .  b2n  y  3/2  »  and 

a2n  b2n+2 
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a2n+2  _  \  2  for  all  sufficiently 

b2n+2  2  b2n 

large  n’s.  If  N  is  one  of  these  n*s,  then  repeated  applications 
of  the  above  inquality  shows  that 

a2N+l 


■i+1  \  (2.)  .  _2N  \ 

/  \2I  h  / 
■J+l  x  x  J  2N 


2h'+l 

for  all  sufficiently  large  k's.  Therefore 

a2n  \  1  ,  or  a„  >  b 

- -  7  2n  7 


2n 


2n 


if  n  is  large  enough.  Thus  y.(X)  does  not  approach  zero  as 


ix2 


4x‘ 


)x|->oo  .  Q.E.D. 

By  a  similar  argument  it  can  be  shown  that  y2(x)e~*A - >0  as 

J  x  |  ~^>cd  |  with  parameters  p=l»3i5»**»  •  Therefore  our  desired 


solutions  are  constant  multiples  of  Hermite  polynomials  HQ(x), 

H^Cx),  H2(x),  ...  .  It  has  been  shown  that  the  Schroedinger 

2  2 

equation  can  be  put  in  the  form  d  z  +  (2p  +  1  -x  )z=0  , 

dx 

which  can  be  transformed  into  a  Hermite  equation  that  has  solutions 


satisfying  2 

ye2*  — ag  | X(  oo  . 

Hence  the  Schroeding  equation  has  solutions  satisfying  w-»0 

as  |u|  cd  if  and  only  if  2E  =  2p+l  or  E=hv(n+£)  where  n  is  non- 

vh 

negative.  Therefore  we  conclude  that  solutions  of 

=  0 


dfw  +  / 2E  -  u2  w  = 

du2  Vhv  / 
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i.„2 


have  solutions  of  the  form  w  -  ce~-’U  H  (u)  where  c  is  a  constant. 

iition  |w|  'du  =  2tfgvnf’  must  be  examined. 

A-  ^  vT 


”aw  the  second  condi 


-l-u2 


CO 


f  c.  _ 

Substituting  w  =  ce  2“  Hn(u)  we  have  J  c2e*"u  ( u )j  2du. 

This  integral  can  be  solved  using  the  orthogonality  property 
when  n =m  and  an  application  of  Roaigues'  formula.  This  formula 
for  the  Kermite  oolynomials  states 

2  2 

Hn(x)  =  (-1)  e  d —  e 


n 


Therefore 


J  c2e“u  [Hn(u)j  2du  =  c2  f  e“u 


dx 


e  Hn(u)Hn(u)  du  . 


Applying  Rodrigues'  formula,  we  have 


°°  2  ( 
e“U  Hn(u)Hn(u)du  =  c2(-l)n 

r*  J 

=  c2(-l)n  I'  c£_  e'^H^ u)du  . 
Integrating  by  parts,  we  have 


•  o o 


J 


-  O 


-  oo 


e"u  eu  d^_  e-u^H  (u)du 
dxn 


J  dx 

-OO 


w=Hn(u) 
dw=Hn' (u) 


,  dv=dn  e-u  du  ,  and 
d^I 


v=d 


n-l  e 


-u 


du 


n-l 


2  2 

Note  that  wv=dn”1e~u  Hn(u)  is  the  product  of  e“u  and  a  poly¬ 
nomials  and  vanishes  at  both  limits.  Therefore 
,o°  'G* 


j e'“2  Oh'"*} 

-  oo 


2du  =  (-l)n+1 


L  /^'(uld—1  e“u 

J L  du-1 


du 


i 
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(_i)n+2  \h  •• 
v  i  n  „  o 

du  2 


( u )dr‘~2  e-u  du 


=  (-D 


2n 


jV" 


^  (u)e  u  'du 


The  term  containing  the  highest  power  of  u  is  2nun»  therefore 

Hn^u^(u)  =  2nn!  and  the  last  integral  becomes 

2  fte  2 

2nn!  )  e“u  du  =  2n  2  n!  (  e_u  du 
Oa 


<3b 


Since 


1  e~u  du  is  a  gamma  function,  we  have 

-'O  ^ 

/  <2°  o 


2nn!  e"u  du  =  (2nn!)2j*r  =  2nn!  . 

2 


Therefore 


-Oo 


H'u2rHn<',>32du  -  2nn!fi?. 


and  the  second  condition  becomes 
o o 


C  r - “ - ] 

I  jwj2du  =  2B\[vtti  =  try  4vm  ,  and 

I  '  h  '  h 

c2  J e-u  [Hn(u)J2du  =  TT\[4vm'1  , 


4  vra 


22n(n!)2h 
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There? ore  the  Schroeding  equation  d_2w  +  h E  -  u2\w  -  0  v/i  th 


du‘ 

,(C 


nv 


(,c  ? 

conditions  w  >0  as  l>:\  — >  go  and  I  Ivn  dx  =  1  is  of  the  form  of  a 
Sturm-Liouville  equation.  The  nontrivial  solutions  exist  if  and 
only  if  the  eigenvalue  X  is  X  ~E=hv(n+J)  where  nlO.  The  non¬ 
trivial  solutions  corresponding  to  the  values  of  are  the  eigen¬ 
functions 


w  = 


4  vm 


22n(n* )2h 


A  _l-u2 
4e  H 


n(u) 


The  concept  of  quanitatively  analysing  differential  equations 
can  be  applied  to  understanding  the  movement  of  a  swinging  pen¬ 
dulum.  For  a  simple  pendulum  one  measures  the  position  of  the 
swinging  arm  by  an  angular  coordinate  0  .  We  can  learn  a  good 
deal  about  the  motion  of  a  pendulum  if  we  convert  the  equation  of 
motion  into  a  system  and  discuss  its  phase  portrait  and  stability 
of  some  of  its  solutions.  It  is  well  known  that  energy  is 
dissipated  in  the  action  of  any  real  dynamical  system,  usually 
through  some  form  of  friction.  However,  in  certain  situations 
this  dissipation  is  so  slow  that  it  can  be  neglected  over  relatively 
short  periods  of  time.  In  such  cases,  we  assume  the  law  of  con¬ 
servation  of  energy,  namely,  the  sum  of  the  kinetic  and  poten¬ 
tial  energy  to  be  constant.  A  system  of  this  kind  is  said  to  be 
conservative.  The  simplest  conservative  system  consists  of  a  mass 
in  attached  to  a  spring  and  moving  in  a  straight  line  through  a 
vacuum.  If  we  let  x  represent  the  displacemnet  of  m  from  its 
equilibrium  position,  and  let  the  restoring  force  exerted  on  ra  ’  ■ 
by  the  string  be  -kx  where  k  >0,  then  the  equation  of  motion  is 
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0 


(3D 


m  d  jk  +  kx  = 

ji.2 


\  prirv:  of  this  >’nd  i  s  called  a  linear  spring,  because  the  ro:f:or- 
i  ng  force  Is  a  linear  function  of  x.  If  m  moves  through  a  resist¬ 
ing  medium,  and  the  resistance  (or  damping  force)  exerted  on  m  is 

-c(dx/dt)  where  c^O,  then  the  equation  of  motion  of  this  non- 

2 

conservative  system  is  m  d  x  +  c  dx  +  kx  =  0  .  Here  we  have 

dt2  dt 

linear  damping  because  the  damning  force  is  a  linear  function  of 

dx  .  If  f  and  g  are  arbitrary  functions  with  the  property  that 
dt 

T 

f ( 0 ) =0  and  g(0)=0,  then  the  general  equation 


m  d*Sc  +  g/dx\  +  f(x)  =  0 
dt2 


can  be  interpreted  as  the  equation  of  motion  of  a  mass  m  under  the 
action  of  a  restoring  force  -f(x)  and  a  damping  force  -g(dx/dt). 

In  general  these  forces  are  nonlinear  and  equation  (31)  can  be 
regarded  as  the  basic  equation  of  nonlinear  mechanics.  In  our 
discussion  we  shall  consider  the  special  case  of  a  nonlinear 
conservation  system  described  by  the  equation 

m  dfx  +  f(x)  =  0  , 
dt2 

Where  the  damping  force  is  zero  and  there  is  consequently  no 
dissipation  of  energy.  Equation  (32)  is  equivalent  to  the  non¬ 
linear  autonomous  system 

dx/dt  =  y  ,  dy/dt  =  -  f(x)/ra 
Eliminating  dt,  we  obtain  the  differential  equation  of  the  paths 
defined  by  solutions  of  (33)  in  the  x,y  phase  plane, 
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'  dy/dx  -  -f(x)/m y  . 

Separating  the  variables  in  equation  (3*0  we  obtain 

mydy  -  -f(x)dx  . 


Su noose  >: -x 
-  *  o 

from  t  to  t 
o 


and  y-yQ  when  t=tQl  then  by  integrating  equation(35) 
we  obtain 


.  my  -  ;;myr 


f (x)ax, 


or 


^x 


r 


+  ]f(x)dx  =  Jmy Q  +  jf(x)dx  . 


2  2 

Now  I  my  =  lra(dx/dt)  is  the  kinetic  energy  of  the  system  and 


v(x)  =  - 


x 

f (x )dx 


<3'0 


(35) 


(36) 


is  the  potential  energy.  Thus  equation  (36)  takes  the  form  of 

Jrmy2  +  v(x)  =  E  ,  (37) 

where  E=;-my0  +  v(x0)  is  the  constant  total  energy  of  the  system. 

Since  equation  (36)  was  obtained  by  integrating  equation  (35) » 
we  see  that  equation  (37)  gives  the  family  of  paths  in  the  x,y 
phase  plane.  For  a  given  value  of  E,  the  path  given  by  equation 
(37)  is  a  curve  of  constant  energy  in  this  plane.  In  other  words, 
along  a  particular  path  the  total  energy  of  the  system  is  a  con¬ 
stant;  this  expresses  the  lav;  of  conservation  of  energy.  The 
critical  points  of  the  autonomous  system  (33)  are  the  points  with 
coordinates  (x  ,0),  where  x  are  the  roots  of  the  equation  f(x)=0. 

C  v 

These  are  equilibrium  points  of  tte  dynamical  system.  From  the 
differential  equation  (35)  we  see  that  the  paths  cross  the 
x-axis  at  right  angles  and  have  horizontal  tangents  along  the  lines 
x=xc.  Equation  (37)  also  tells  us  that  the  paths  are  symmetrical 
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with  r  or  poet  to  the  a -ax  '  s.  Solving  coua  ■  '  on  (j7)  for  y,  ,ve  find 

y  ~  +..'2  p>v(x)] 

Vm 

'low  let  us  consider  a  simple  pendulum  which  is  composed  of  a 
mass  m  at  the  end  of  a  straight  wire  of  negligible  mass  and 
length  n.  fuorose  this  pr  ndulum  is  s.u  ncnled  from  a  f i y ed  point 

and  is  free  to  vibrated  in  a  vertical  plane.  Let  x  denote  the  angle 
which  the  straight  wire  makes  with  the  vertical  ray  down¬ 
ward  at  time  t.  Assume  there  is  no  air  resistance  and  the  only 
forces  acting  on  the  mass  are  the  tension  in  the  wire  and  the 
force  due  to  gravity.  The  differential  equation  for  the  displace¬ 
ment  of  the  pendulum  is 

2 

d  x  +  ksin  x  =  0  , 
dt2 

where  k  is  a  positive  constant.  The  above  equation  is  in  the 
form  of  equation  (32),  where  m=l  and  f(x)=ksin  x.  The  autonomous 
system  equivalent  to  equation  (38)  is 

dx/dt  =  y  ,  dy/dt  =  -  ksin  x  . 

The  critical  points  of  this  system  are  points  (x  ,0),  where  x 

are  the  roots  of  sin  x=0.  Thus  the  critical  points  are  an 

infinite  set  of  points  (nTT^O)  where  n=0,+l,+2, . . .  .  Hence  the 

differential  equation  of  the  paths  in  the  x,y  phase  plane  is 

dy  sin  x  . 

dx  y 

By  variable  separation  we  have 


(38) 


ydy  =  -  sin  xdx 


?o  :.he 


equation  of  the  family  of  paths  is 


hy2  +  (k-kcos  x)  =  E  . 

2 

'•'qua '.ion  (.19)  5  s  of  the  form  of  equation  (37)  where  ia=l ,  ;’.y 

is  the  kinetic  energy,  and  v(x)  =  f(x)dx  =  k-kcos  x  is  the 

7o 

potential  nnerry.  Mow  we  shall  construct  the  xy  phase  plane 
dia  ram  by  first  constructing  the  curve  z~v(x)  and  the  lines 
z=E  for  various  values  E>0  (see  Figure  1,  where  z=E=2k).  From 
this  plane  the  values  E-v(x)  can  be  read  on  the  graphs  and  can  be 
plotted  in  the  xy  phase  plane  directly  below  the  value  of  y  as 

y  =  + 

=  +  v^2  ^E-(k-kcos  xjj 

Using  the  formula  v(x)=k-kcos  x  for  the  potential  energy,  then 

v'(x)=  sin  x  and  v*’(x)=  cos  x  . 

However,  v*(x)=0  at  x=n7T  where  n=0,+l,+2,  ...  .  Observe  that 

v,'(nIT)>0  if  n  is  even  and  v,'(nTT,)<0  if  n  is  odd.  Therefore 
the  potential  energy  v  has  a  relative  minimum  at  x=2nf>  for  every 
even  integer  2n  (where  n=0,+l,+2,  ...  )  and  a  relative  maximum 

at  x=(2n+ljft  for  every  odd  interger  (2n+l),  where  n=0,+l,+2,  ...  . 

Hence  the  critical  points  (2nfT,0)  where  n=0,+l,+2, . . . )  are  centers 
and  are  stable;  and  critical  points  ^2n+l),oJ  fahere  n=0,+l,+2,  ...) 
are  saddle  points  and  unstable. 

(Note;  The  diameter  of  a  closed  curve  c  is  called  the  maxi¬ 
mum  of  the  distance  (A,B)  between  the  points  A  and  B  on  the 
curve  c  for  all  possible  pairs  of  points  A  and  B  on  c.  Con¬ 
sider  Figure  3»  the  critical  point  (0,0)  is  called  a  saddle 
point  and  such  a  point  may  be  characterized  as  follows; 
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(39) 


1.  It  is  approached  and  entered  by  two  half-line  par/ns 

(AO  and  -50 )  as  t  ">oo  .  These  two  paths  forming  the  roome trie 
curve  AB. 

2.  It  is  approached  and  entered  by  two  half-line  paths 
(CO  and  DO)  as  t.  xd  .  These  two  paths  forming  the  gc-o- 
inetric  curve  CD. 

3.  Between  the  four  half-line  paths  described  in  (1)  and 

(2)  there  are  four  domains  R^,  Rg,  R^,  R^,  each  containing 
an  infinite  family  of  semi-hyperbolic-like  paths  which  do 
not  approach  0  as  t->+  oo  adn  t“*-oo  but  which  become  asympto¬ 
tic  to  one  ■  of  the  four  half-line  paths  as  t-^+co 

and  as  t-^v-ao  .  ) 

Using  Figure  2  we  can  see  that  if  the  total  energy  E  is  less 
than  2,  then  the  corresponding  paths  are  closed.  Each  of  these 
paths  surrounds  one  of  the  centers  (2nlt,0).  Physically  all  of 
of  these  stable  centers  correspond  to  exactly  one  physical  state, 
namely,  the  pendulum  at  rest  with  the  bob  in  its  lowest  position 
(the  stable  equilibrium  position).  Therefore  each  closed  path 
about  a  center  corresponds  to  a  periodic  back-and-forth  oscil¬ 
lation  motion  of  the  pendulum  bob  about  its  lowest  position. 

On  the  other  hand,  if  the  total  energy  is  greater  than  2,  then  the 
corresponding  paths  are  not  closed.  Clearly  x->+cd  as  t«^+a> 

if  y-dx >  0  ,  and  x— >-oo  as  t-*-oo  if  y=dx<  0.  Thus  the  motion 
dy  dy 

corresponding  to  such  a  path  does  not  define  x  as  a  periodic 

function  of  t.  Nevertheless,  physically,  the  corresponding  motion 

is  periodic.  To  see  this  observe  the  following!  All  the  unstable 
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corners  (( ?n+l ) , o)  corespond  to  exactly  one  physical  state, 
namely,  the  pendulum  at  rest  with  the  bob  in  its  highest  position; 
and  each  of  the  nonclosed  paths  corresponds  to  a  physical  rota¬ 
ting  motion  of  the  pendulum  about  its  pivot  point.  If  the  total 
on -Try  value  is  '1--2 ,  this  separates  the  bad:  -  and  -  for  th  motion  from 
the  rotational  motion.  The  paths  that  correspond  to  are  the 

ones  which  enter  the  saddle  points  as  t -^ao  or  t  =>  -a>  . 

Now  we  turn  our  discussion  to  oscillations  of  second  order 
nonlinear  equations.  Recall  an  n^-order  ordinary  equation  in 
the  unknown  y  and  independent  variable  x  is  nonlinear  if  it  can 
not  be  put  in  the  form 

t>n(x)d^y  +  bn  1(x)dn~1  y  +  ...  +  b^xjdy;  +  bQ(x)y  =  g(x). 
dxn  dx11-1  dx 

The  well  known  Van  der  Pol  equation  is  a  second  order  nonlinear 

equation.  This  equation  is 

d2x  +  E(x2-1 )dx  +  x  =  0  , 
dt2  dt 

where  E  is  a  positive  constant.  We  shall  investigate  the  Vain 
der  Pol  equation  and  see  if  there  exist  periodic  solutions. 

Theorem  ?  can  be  used  to  determine  if  the  Van  der  Pol  equation 
has  periodic  solutions  providing  all  of  the  requirements  can  be 
met.  Checking  the  properties  of  the  theorem,  we  have 

1.  Since  f(-x)=E(x2-l)=f(x) ,  then  the  function  f  is  even. 

Clearly  f  is  continuous  for  all  x. 

2.  For  F(x),  we  have  F(x)=E(x3-x) .  Note  if  x=  /S',  then 

3 

F(x)=EQ3/2-3*)  =  E(33//2-33//2)  =  E(0)  =  0 

3  3 


ko 


pherefore  for  04  x  <  {3  ,  then  F(x)-E(x^-x)  is  negative. 

,  -  3 

For  x^,1!,  then  F(x)  is  positive.  In  fact,  F(x)  is  mono¬ 
tonic  increasing  and  F(x)  >co  as  x  -*oo  . 

J.  Since  g(-x)  =  -x  =  -g(x),  if  follows  that  the  function 
r  is  old.  Also  since  g'(x)--l,  then  the  derivative  of  g  is 
continuous  for  all  x.  It  is  obviously  true  that  g(x)^>  0 
for  x  >0. 

Clearly  G(x)->ao  as  x^>co  . 

'nhus  the  Van  der  Pol  equation  satisfies  all  the  properties 
of  the  Lienard-Levinson-Smith  Theorem  and  therefore  has  a  unique 
nontrivial  periodic  solution.  In  other  words,  the  equivalent 
autonomous  system 

dx/dt  =  y  ,  dy/dt  =  E(l-x2)y  -  x 

has  a  unique  closed  path.  Reducing  the  above  autonomous  system 

to  a  differential  equation  of  the  paths  we  have 

dy  dx  =  (i:(l-x2)y  -  xl  •  y 
dt  dt 

Using  the  chain  rule 

dy  dt  =  fe(l-x2)y  -  x]  1 
dt  dx  y 

dy  =  E(l-x2)y  -  x 
dx  y 

Graphs  of  the  Van  derPol  equation  appear  differently  than 
graphs  of  equations  that  satisfy  the  Stum  Comparison  Theorem. 
However,  they  depict  essentially  the  same  meaning.  A  graph  of 
an  equation  that  satisfy  the  Sturm  Comparison  Theorem  appear  as 
in  the  following  figure. 
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* 


The  equation  oscillates  about  the  x-axis  and  reaches  its 
limit  about  that  axis.  Graphs  of  the  Van  der  Pol  equation  approach 
a  limit  cycle.  This  cycle  is  unique  whether  the  approach  is  made 
from  inside  the  closed  curve  or  outside.  Graphs  of  equations  that 
satisfy  the  Sturm  Comparison  Theorem  will  resemble  graphs  of  the 
Van  der  Pol  equations  if  one  would  stand  at  the  end  of  the 
x-axis  and  look  down  the  x-axis  toward  the  origin.  In  this 
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perspective,  one  could  observe  the  path  taken  to  reach  the  limit 
cycle . 

nhe  following  are  graphs  of  the  Van  der  Fol  equation  with 
different  values  of  E.  The  closer  the  value  of  E  is  to  zero 
the  more  the  graph  approximates  a  circle.  Figure  1  is  a  graph 
with  E-.002.  Even  though  this  graph  is  not  a  perfect  circle, 
it  is  closer  than  any  of  the  others.  Figure  2  shows  how  the 
limit  cycle  is  reached  from  inside  and  outside  the  closed  curve 
with  E=.25  •  The  further  away  the  value  of  E  is  from  zero  the 
more  the  closed  curve  becomes  oblonged  about  the  y-axis.  Figure 
3  also  shows  the  limit  cycle  from  inside  and  outside  the  closed 
curve  with  E=1.5*  Figure  4  shows  the  limit  cycle  with  E=4.  Note 
how  the  graph  is  becoming  more  oblonged  and  the  values  of  y  are 
increasing.  Figure  5  shows  a  more  oblonged  graph  with  E=7* 

Figure  6  shows  the  differences  in  the  limit  cycles  with  all  the 
graphs  starting  from  the  same  point,  x=1.5i  and  different  values 
of  E. 
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